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Abstract 

We study the cl = 25 limit, which corresponds to c = 1 string theory, of bulk and 
boundary correlation functions of Liouville theory with FZZT boundary conditions. This limit 
is singular and requires a renormalization of vertex operators. We formulate a regularization 
procedure which allows to extract finite physical results. A particular attention is paid to 
c = 1 string theory compactified at the self-dual radius R = 1. In this case, the boundary 
correlation functions diverge even after the multiplicative renormalization. We show that all 
infinite contributions can be interpreted as contact terms arising from degenerate world sheet 
configurations. After their subtraction, one gets a well defined set of correlation functions. We 
also obtain several new results for correlation functions in Liouville theory at generic central 
charge. 
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1 Introduction 

Non-critical string theories (see for instance [1-3] for a review) have long played the role of an 
interesting laboratory to test various phenomena of critical string theories. A recent breakthrough 
in the understanding of their non-perturbative effects [4-23] drew a lot of attention to these theories. 
The ground for this breakthrough was prepared by a progress in Liouville theory, where two sets 
of conformally invariant boundary conditions, corresponding to the so called FZZT [24,25] and ZZ 
branes [26], were discovered and all basic boundary correlation functions were calculated [27-29]. 
A review of these results can be found in [30]. 
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Usually, the Liouville correlations functions are found for generic values of the Liouville central 
charge and they are represented as complicated combinations of (integrals of) some special func- 
tions. On the other hand, the most interesting non-critical string theory is the c = 1 string theory. 
Besides giving the first non-trivial example of a non-rational CFT, it possesses a clear geometric 
interpretation as a string theory in a two-dimensional spacetime. It also allows for non-trivial 
backgrounds, such as the two-dimensional black hole, so that it might be useful in the analysis of 
black hole thermodynamics. Finally, it appears as a building block in many constructions of critical 
superstrings. 

The field content of c = 1 string theory consists of a massless scalar matter coupled with the 
Liouville field with central charge cl = 25. This is a very special value of the central charge for 
Liouville theory, since its correlation functions typically diverge at this point. This well known 
divergence is a counterpart of the fact that the c = 1 CFT is a singular limit in the family of 
minimal (p,q) models [31]. 

In this paper we demonstrate how it is possible to define all Liouville correlation functions in 
the Cl = 25 limit systematically. The above mentioned divergence at cl = 25 can be removed by a 
multiplicative renormalization of the Liouville coupling constants. Whereas this is a straightforward 
procedure for the theory in the bulk, for the boundary correlation functions there are additional 
complications related to the appearance of the so called "contact terms" , which diverge even after 
the renormalization. In the simplest cases, these contact terms are entire functions of the boundary 
cosmological constants and therefore vanish after the insertion of a finite number of boundary 
cosmological constant operators. In the matrix model approach, such terms are usually considered 
as non-universal and excluded from the final results, and one can expect the same situation to 
occur in the CFT under consideration. 

However, our analysis shows that there are other divergent contributions, which have a more 
complicated form and have no simple analogues in the matrix models. Nevertheless, we show 
that all of the divergent terms arising in the cl = 25 limit, independently of their form, can be 
interpreted as contributions coming from degenerate two-dimensional geometries, and they should 
therefore be excluded from the resulting expressions for the correlation functions. In this way one 
arrives at finite results for all basic correlators, which, once multiplied by the corresponding matter 
contributions, give the correlation functions of c = 1 string theory. 

Recently, an attempt to extract the cl = 25 limit of Liouville theory with boundary was made 
in [32] . However, the complication related to the appearance of the contact terms was not taken into 
account. As a result, [32] provided only divergent results for the leading non-universal contributions, 
whereas the physical part of the correlation functions remained non calculated. Our paper fulfills 
this gap and presents the physical results derived from a careful regularization procedure based 
on a consistent coupling of Liouville theory and matter CFT. As in [32], we limit our analysis to 
boundary conditions of FZZT type. The correlation functions on ZZ branes can be obtained by 
applying a well known relation between these two types of branes [5, 11, 27, 29, 33]. 1 

The singularity of the Cl = 25 limit implies that the final results are simpler than those for 
generic values of the central charge. The greatest simplifications occur if one restricts oneself 
to the spectrum of vertex operators corresponding to c = 1 string theory compactified at the 
self-dual radius R = 1, which is the case considered in [32]. The reason is that in this case 
all operators are degenerate operators in both the matter and Liouville sectors. In particular, 
all special functions reduce to ordinary functions and all correlators can be given quite explicit 

: Let us note that since the results of [32] depend analytically on the boundary cosmological constant, this relation 
implies that the corresponding correlation functions with ZZ boundary condition vanish, which is clearly unphysical. 
In contrast, our results will ensure non-vanishing correlators for ZZ branes. 
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expressions. Therefore, we will be primarily interested in this situation. The self-dual radius 
introduces additional singularities with respect to the generic case. We show that, after subtraction 
of the contact terms, our renormalization ensures that all correlation functions remain finite. 

The restriction to the self-dual radius of compactification is also relevant for the link to topolog- 
ical strings. It has long been known that c = 1 non-critical string theory compactified on a circle of 
radius R = 1 provides a dual description of topological strings on the conifold [34-38]. One might 
therefore expect that the correlation functions that we study here shed light on the dynamics of 
some D-branes in topological string theory. 

The organization of the paper is the following. In section El we recall some basic facts about 
Liouville theory, introduce our notations and collect the known results for the correlation functions 
of Liouville vertex operators on a sphere and on a disk. In section |3] we formulate our procedure 
to take the c = 1 limit, based on a CFT representing Liouville theory coupled to matter with 
background charge. At this point we also introduce the multiplicative renormalization of the 
vertex operators. In the following two sections we compute sphere and disk correlation functions 
for Cl = 25. The former are considered in section 4 and the latter are dealt with in section 
5. Since in the case of boundary correlators the multiplicative renormalization is not enough to 
remove all divergences, we begin in subsection 15.11 by discussing the meaning of the contact terms 
and the necessity of excluding them from the results. We then proceed with the details of the 
computation in the remaining part of the section. Notice that, as an intermediate step in our 
derivation of the correlators in the cl = 25 limit, we also find the Liouville correlation functions for 
a generic central charge with momenta corresponding to degenerate matter fields, and these agree 
with the predictions of the microscopic approach found in [39]. Finally, we conclude in sectional 
by summarizing the list of our main results. Readers who are not interested in the details of the 
derivations may find this list useful to get to the results of our work quickly. Several appendices 
present various details of our calculations. 



2 Correlation functions in Liouville theory 

Liouville theory arises when one considers two-dimensional quantum gravity, possibly coupled to 
matter, in the conformal gauge. It is defined by the following action: 

It is a CFT whose central charge is given by 

c L = l + 6Q 2 (2.2) 
and the parameter b in (j2.1j) is related to Q via the relation 

Q = b+\. (2.3) 
b 

In general, b and Q are determined by the requirement that the total central charge of the Liouville 
and matter field is equal to 26. For example, in the case when matter is represented by a minimal 
(p, q) model with central charge c Pi(J = 1 — 6 ^~^ , the relation ()2.2|) implies that b = \Jp/q. 
Coupling to c = 1 matter corresponds to the limit b = 1 so that the Liouville central charge 
becomes cl = 25. 
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2.1 Sphere correlation functions 

An important class of conformal primaries in Liouville theory consists of the bulk operators 

V a (z) = e 2a ^ z) , (2.4) 

whose scaling dimension is given by A a = A a = a(Q — a). In particular, the Liouville interaction 
in the action ()2.1|) can be seen as a perturbation of the Lagrangian by the cosmological constant 
operator SC = [iV^. 

In principle, all information about Liouville theory on a manifold with topology of a sphere is 
contained in the three-point correlation function, which was found to be [40,41] 

. Ct 2 , C3) 

(V ai (zx)V a2 {z 2 )V a3 (£ 3 ))sphere = j^j 2 (A 1+ A 2 -A 3 ) | 2(A 2 +A 3 -Ai) ^ | 2(A3+A!-A 2 ) ' ( 2 " 5 ) 

where 

C(a 1 ,a 2 ,a 3 )= [nw{b 2 )b 2 - 2b \ (2.6) 

x T b , T b (2a 1 )T b (2a 2 )T fe (2a 3 ) 

Tfi(ai + a 2 + a 3 — Q)T b (a 1 + a 2 — a 3 ) Y 6 (a 2 + a 3 — aa) Y b (a 3 + ai — a 2 ) 

Here we used the standard notation for the di-gamma function 

**> = rKj- (2 ' 7) 

while the other special functions appearing in ()2.6)1 and in the following are defined in appendix lAl 
A particular important case of the three-point function is given by the two-point function. In 
Liouville theory it is also often interpreted as a reflection amplitude which is the coefficient Z)( re fl)(a) 
appearing after the duality transformation a — > Q — a: 

V a (z) = D (lc& )((x)V Q - a (z). (2.8) 

Its explicit expression reads 

n f \ r ^2xi(Q-2a)/5 7(26a-6 2 ) 

However, in string theory and quantum gravity the correct way to define the two-point function 
is through the three-point correlator (|2.6|) . As explained for instance in appendix C of [42] 2 one 
should take the three-point function at the particular values of momenta a.\ = a 2 and a 3 = b, and 
then integrate it once with respect to — \x. This procedure allows to avoid any problems related to 
the residual gauge symmetry on the sphere with two insertions of the vertex operators. In this way 
one obtains 

r>/ x r^.„.,^2M(Q~2a)/b (g - 2a)7(25a - b 2 ) 

D(a) = [ W (6 )] 7r6 2 7(2 _ 2a/6+1/6 2 ) - ( 2 - 1Q ) 

Thus, the two-point function found from the three-point correlator differs from the reflection am- 
plitude by a factor: 

D{a) = 9_^ D {a) . (2.11) 

In our study we have to use the function (|2.10j) . In particular, after the renormalization only D(a) 
gives finite expressions for all momenta in the c = 1 limit. 

2 We thank Joerg Teschner for drawing our attention to this reference. 
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2.2 Disk correlation functions 



When Liouville theory is formulated on a surface with boundary, there are two possible boundary 
conditions, respectively of Neumann and Dirichlet type. The first type gives rise to the so called 
FZZT branes [24,25], while the second one leads to the ZZ branes [26]. In order to work with FZZT 
boundary conditions, one adds the following boundary contribution to the action (|2.1j) : 

Shnd = 27 J 9E ( Qk<P + 27rfi ^) ■ (212) 

The boundary conditions are parameterized by the boundary cosmological constant fi B . However, 
all correlation functions are more naturally written in terms of another parameter s, which is related 
to fi B through 

a*b = J—n^ cos H* bs )- ( 2 - 13 ) 

V sm(7rtr) 

In this theory, besides the bulk operators (j2.4j) . there are operators acting on the boundary and 
changing the boundary conditions 

= e^tt>. (2.14) 

Here the labels si and s 2 correspond to the boundary conditions to the left and right of the point 
£ on the boundary. The scaling dimension of the operator ()2.14j) is given by A^ = A^ = [3(Q — l3). 

Let us now summarize the non-trivial correlation functions in the theory with boundary. We 
recall that definitions and properties of the special functions appearing below are summarized in 
appendix |XJ 

• One-point bulk correlation function [24] 

(V a (z)) disk = UH " ) A , (2.15) 



U(a\s) = \ YR^{b 2 )] iQ - 2a ^ 2h Y{2ba - b 2 )T(2a/b - 1/b 2 - 1) cosh ((2a - Q)irs) . (2.16) 
Bulk-boundary correlation function [27] 

(v a (z) B -(ou sk = lz _- l2 fX\ s L^r^ (2 - 17) 



R(a,[3\s) = -2m n^(b 2 )b 



2-2*1 ^- 2a -^ 2b T 3 b (Q - (3)T b (2a - (3)T b (2Q - 2a - l3) 



dt e 



-2nts 



T b (Q)T b (Q - 2f3)T b ((3)T b (Q - 2a)T b (2a) 
S b (t + (3/2 + a- Q/2)S b (t + (3/2-a + Q/2) 



S b (t -P/2-a + 3Q/2)S b (t -f3/2 + a + Q/2) ' 



(2.18) 



• Three-point boundary correlation function [28] 

(Bp l 3 S3 (Ci)Bg S2 ^ 2 )B^ s 1 (6))disk = , , Ai + A ^f 3 1 f | if+ i^l' | ? n A3+A1-A2 ' ( 2 - 19 ) 
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C(ft,ft,ft|si,s 2 ,s 3 ) = -Am [7r/i 7 (&> 2 - 2b 

x T b ((3 2 + /3 3 - ft)r fe (2Q - ft - ft - ft)r b (Q - ft - ft + j3 3 )T b (Q - ft + ft - ft) 

T b (Q)T b (Q - 2ft)r b (Q - 2(3 2 )T b (Q - 2ft) 



S 6 (Q-ft + i^)S 6 (Q-ft-i^) ( AS^W) 



X 



S 6 (ft + z^a)S 6 (ft-z^±a) 7 fi S *>(* + ^) 



where 



Ul 


= Q + i s -^- 


ft, 




= 2Q + *^-ft-ft, 


U-2 


= Q-t^- 


ft, 


r 2 


= g + ^-ft + ft, 




= z^+ft, 




V3 


= Q + is 2 , 


u 4 


= Q + i 8 -^- 


ft, 


^4 


= Q. 



(2.21) 



Several comments are in order. First of all, one should give a precise definition of the contours 
of integration in (|2.18|) and (|2.2U|) . The rule is the following: for both integrals, the integration 
contour in the complex t-plane lies on the right of the poles of the integrand coming from poles 
of the numerator and to the left of the poles coming from zeros of the denominator. This gives 
an unambiguous definition of the integrals, except in the case when a pole from the numerator 
coincides with a pole from the denominator. If such a situation of colliding poles occurs, it means 
that the integral has a singularity at the corresponding point in the parameter space. 

Another comment is that we have changed the overall normalization of the one-point bulk 
and the three-point boundary correlation functions with respect to the one usually found in the 
literature. Compared to the standard expressions from [24] and [28], the first one is divided by 2n 
and the second one is multiplied by Air. This is done because the theory is completely defined by 
the bulk-boundary and the three-point correlation functions, whereas the one-point bulk correlator 
is a derived quantity, whose normalization can be fixed starting from any of the two. For example, 
one can take the bulk-boundary correlation function ()2.18|) with the boundary operator being the 
boundary cosmological constant operator, i.e. (3 = b, and then integrate with respect to — /x s . In 
appendix [BJ we show that this procedure gives precisely (|2.16|) with the factor l/2n included. This 
factor can be traced back to the rotational symmetry of the disk, which remains non-fixed in the 
one-point correlator. The correct normalization valid for string theory requires to divide by the 
volume of the gauge group and thus follows from the bulk-boundary correlation function where all 
symmetries are fixed. 

However, a similar procedure (see below) done for the standard expression for the three-point 
function leads to a result differing by a factor An. This means that there is a mismatch between 
the bulk-boundary and the three-point functions. Assuming that the bulk-boundary correlator is 
correct, one should then modify the standard expression for the three-point correlator by including 
this factor, as we did in (|2.2U|) . 

The third comment is that, similarly to the bulk case, one can also consider the two-point 
boundary correlation function. Its usual definition used in Liouville theory is through the reflection 
relation [24] 

S? -a (0=d(»a)C9|ai,a2)B55(0, (2.22) 

where 

r ,i (Q-2/3)/2b 

\K m {b 2 )b 2 - 2b2 \ T b {2f3 - Q)T-\Q - 2(3) 

d(refi)(ftsi, s 2 ) = Sfe ^ + 2 £i±£2) Sb _ pT±^) Sb ({3 + i^ffi) S b {(3 - i 3 -^)' (2 ' 23) 
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However, as usual, the two-point function relevant for quantum gravity should be obtained through 
the three-point function (|2.20|) by the standard procedure of taking one of the operators to be the 
boundary cosmological constant operator and integrating with respect to —\i B . In fact, since the 
initial expression is not explicitly symmetric, there are three inequivalent ways to do this depending 
on which operator we choose to represent B^ s . In appendix O we show that all three ways lead 
to the same result. This fact provides a non-trivial explicit check of the symmetry of the three- 
point function ()2.20|) under cyclic permutations of its arguments. We find that, as in the bulk case 
(|2.1ip . the two-point function differs from the reflection amplitude ()2.23|) by a simple (momentum 
dependent) factor 

d(P\s 1 ,s 2 ) = 2(Q-2P)d (re{l) (P\s 1 ,s 2 ). (2.24) 

In appendix 10.31 starting from the two-point function found in this way, we also rederive the 
one-point function of the cosmological constant operator. The result should coincide with the 
one obtained from the bulk-boundary correlator, and this forced us to introduce the An factor 
mentioned above. These results also contain information about the one-point function of the 
boundary cosmological constant operator which is found to be 

lb i / i\ 

W{s) = — [tt//7(6 2 )] T f 2 - - J cosh(vrs/6). (2.25) 

This quantity will be quite important in the following so it is useful to give its explicit expression 
at this point. 

3 The c = 1 limit 

3.1 2D string theory with a background charge 

It is well known that the c = 1 limit of the family of CFTs with central charge c < 1 coupled 
to Liouville theory with cl = 26 — c is singular. Therefore, one cannot put directly b = 1 in the 
results obtained in Liouville theory with a generic central charge. To make the limit meaningful, 
we perform a regularization by considering Liouville theory coupled to a Gaussian matter field with 
a background charge q, following [39,43]. The Gaussian field is described by the following action 

S mat = y I ' d 2 z^ ((&r) 2 + iqxR\ + ±- ! (iqxK) . (3.1) 

The background charge leads to a shift of the central charge of the matter field 

Cmat = 1 - 6g 2 . (3.2) 
Requiring that the total central charge of the system vanishes 

Ctot = c L + c mat + c ghost = (1 + 6Q 2 ) + (1 - 6q 2 ) - 26 = 0, (3.3) 
one can express the background charge in terms of the Liouville parameter b 

q = \-b. (3.4) 
The bulk and boundary vertex operators of the full theory 

V p = e-^e^W B p = e - ipx e m * (3.5) 
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are found from the condition that their dimensions are equal to 1 

A[V P ] = |(! + ?)+a(Q-a) = l, 
A[B P ] = p(p + q)+p(Q-l3) = l. 

Taking into account the Seiberg bound [44], this gives for the parameters a and f3 



a(p) = b — p/2 if p > —q, a(p) 
P(p) = b-p if p> -q/2, pip) 



l/b + p/2 if p < -q, 
1/b + p if p < -q/2. 



(3.6) 
(3.7) 



{3.1 



The two branches of the solution correspond to the left and right moving (or positive and negative 
chirality) tachyons in the target space picture of two-dimensional string theory. 

In order to restrict ourselves to the spectrum of the theory compactified at radius R, we have 
to take the discrete set of operators with quantized momenta 



Pn 



n/R. 



Then for n > the corresponding Liouville momenta are 



a,. 



h n 
b ~2R' 



Pn 



b ~R- 



(3.9) 



(3.10) 



In fact, for generic R and n ^ 0, the correlation functions do not exhibit singularities. Therefore, 
the details of the limiting procedure are not important in this situation. 

However, if R is rational some additional singularities appear so that the c = 1 limit should be 
treated with a great care. We will pay special attention to the case of the self-dual radius R — 1. 
On the one hand, it is the most singular case, but on the other hand it allows the most explicit 
representation for the correlation functions. The origin of the additional divergences can be traced 
back to the fact that all operators with momenta p n = n are Liouville dressed degenerate matter 
fields of the c = 1 CFT. Therefore, it is natural to take the c = 1 limit in such a way that the 
operators in question are degenerate for all b. For general b the momenta corresponding to the bulk 
degenerate fields are 



bulk 



a, 



Q/2 



—q + r/b — sb, 
and those of the boundary degenerate fields are 

pJ? = (-<? + r/b - sb) /2, p r , s = Q/2 



r sb 
2b ~ ~2 



r sb 
2b ~ ~2 



(3.11) 



(3.12) 



where in both cases r, s G Z, rs > 0. Thus, when dealing with the case R — 1, in the regularized 
theory we will set R = b so that p n = n/b and the Liouville momenta ()3.10|) take the form 



OC n = - — , Pn = b-~, 



n > 0. 



(3.13) 



This corresponds to the choice r = n + 1, s = 1 for the bulk and r = 2n + l, s = 1 for the boundary 
momenta. 

However for n < this solution does not work. To get degenerate fields with momenta of 
the form (|3.9jl . one has to choose R — 1/6 and, for example, in the bulk case r = l, s = 1 — n 
so that p n = nb. This reflects difficulties in compactifying a field with a background charge. 
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After compactification, there are degenerate fields only of either positive or negative chirality. On 
the other hand, to write non-vanishing correlation functions, one needs to satisfy the momentum 
conservation law in the matter sector, which due to the presence of the background charge q acquires 
an additional term. For the sphere and disk correlation functions it takes respectively the following 
form 

k k I 

5>* ulk = -2q, J>™ k + 5> bnd = -q. (3.14) 

i=l i=l i=l 

This implies that at least one momentum must have r, s < and hence must not be of the form 
()3.9j) required by compactification. Thus, strictly speaking, we have a compactified theory only in 
the case 6 = 1. 

Usually, there are several ways to satisfy the conservation condition (j3.14j) which differ by the 
choice of chiralities of the vertex operators. In the 6=1 limit they correspond just to different 
choices of signs of momenta and should lead to the same results. We are then free to choose a simple 
convention consisting in taking all operators except one to be of positive chirality. In addition, on 
the disk we will always choose a boundary operator to be the one with negative chirality. Then, if 
for positive chirality p^uik = m/b and p\ nd = rrii/b, the remaining momentum must be of the form 

pbulk = _ 2g _ p bnd = _ q _ m / bj ( 3 . 15 ) 

where n = n i an d m = Yli=i n « + Si=i m *' respectively on the sphere and on the disk. The 

momenta f)3. 15j) are also degenerate and can be obtained from (|3.11|) and ()3.12|) for r = — (n+1), s = 
— 1 and r = —(2m + 1), s = —1, respectively. Starting from the momenta ()3.15|) for the matter 
part, one arrives at the same Liouville momenta as in 1)3.13)1 . 

We conclude that, after fixing all Liouville momenta to be of the form ()3.13)) . the only remaining 
free parameter in our system is b. We are interested in the limit b — > 1 where we should find the 
compactified c = 1 string theory. In order to take this limit, we put b = 1 — e for a small parameter 
£ and extract all singular and finite terms in the correlation functions. We will explain that the 
singular contributions are given by the so called contact terms, which are to be neglected, while 
the physical answers are only contained in the finite part of the correlation functions. The exact 
meaning of the contact terms will be explained in section 15.11 



3.2 Renormalization of couplings and vertex operators 

It is known that to get a finite sensible c = 1 limit, one should simultaneously renormalize the 
couplings of the theory (see for example [5,45]). The basic Liouville couplings are renormalized as 
follows 

7T jJj 

fi = 7r 7 (6 2 )/i, fi B = — ^ b2 y (3.16) 

The normalization is chosen so that the parameterization of the renormalized boundary cosmological 
constant looks simply as 

fi B = cosh(r), (3-17) 

where we introduced a new parameter 

r = nbs. (3.18) 

The introduction of the parameter r is very convenient. It allows to avoid complications related to 
the fact that, although in the c = 1 limit the parameter b disappears from (j3.18j) . it is important 
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to keep track of it to correctly derive the physical part of the correlation functions in the boundary 
theory. Another useful notation, which we will extensively use, is 



P*j =V B ( T i) -p. B {Tj). (3.19) 

The definitions ()3.16|) imply the renormalization of the bulk and boundary cosmological constant 
operators. We will show that, in order to get finite correlation functions at R = 1, one should 
renormalize also all other vertex operators. Generalizing (13.16(1 and following [39,43,46], where the 
multiplicative renormalization was used to compare the CFT results with results of the discrete 
approach to 2D quantum gravity, we define the renormalized bulk and boundary operators with 
momentum p as follows 

v,= ' ,. 4 . v„ ^(WW) 

where the Liouville momenta a(p) and (3(p) were defined in (13.8(1 . Their explicit expressions depend 
on the sign of p and thus operators of positive and negative chirality have different renormalization 
factors: 



Notice that the renormalization factors for different chiralities become the same in the 6 = 1 limit 
when expressed as functions of Liouville momenta. 

The renormalized correlation functions are defined in terms of vertex operators V p and B p 

a (or k }, {pf d }\{ Sj }) = ^nv(^n^'fe)) ( 3 - 23 ) 

and can be obtained from the non-renormalized expressions by multiplying them with the factors 
in (J3.2()() . The correlation functions of the c = 1 theory are found then in the limit b — > 1 

A c =i ({rf ulk }, d }\{ Sj }) = limi ({?r k }, {^ bnd }l{^}) • ( 3 - 24 ) 

Actually, in this paper we will consider only the Liouville part of the correlation functions. However, 
especially in the boundary case, this is the main non-trivial part which contains all information 
about the dependence on the boundary cosmological constant. 

It is clear that in the limit b — *■ 1 the renormalization factors of both bulk and boundary 
cosmological constants appearing in ((3.16(1 vanish, whereas for p n with n > (see ((3.21(1 ) the 
corresponding factors diverge as £ _1 in the bulk case and remain finite for the couplings of the 
boundary operators. The necessity of a renormalization with such properties can be understood 
as follows. First, the analysis of sphere correlation functions fixes the renormalization of all bulk 
operators, which can also be confirmed by the fact that it gives a finite result for the one-point bulk 
correlator on a disk. Then, the renormalization of the boundary operators follows from the study 
of the bulk-boundary function and can be verified by analyzing the boundary two- and three-point 
functions. In the next sections we will study in detail all basic correlation functions of Liouville 
theory following this sequence. 
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4 Sphere correlation functions in the c = 1 limit 

First, we evaluate the ci = 25 limit of the bulk correlation functions for Liouville theory on a 
sphere. The main quantity to consider is the three-point function ()2.6|) . For the momenta ()3.10|) it 
takes the form 

\ \ n 2\ 7 2-2fe 2 l ^ + 2^fl^ nz ~ 2 

C(a ni ,a n2 ,a n3 ) = 717x7(6 )6 (4.1) 

T b (b)T b (26 - f ) T fe (26 - f ) T, (26 - f ) 



x 



T 6 (26 - I - a±H±M) t 6 (6 + n ^ 1 R n2 ) T b (6 + Y fc (6 + 



This expression can be simplified if one takes into account the momentum conservation ()3.14|) 
from the matter sector, which implies that = ri\ + 712. We also apply here the multiplicative 
renormalization ()3.2()j) . As explained in the previous section, our convention is to choose two 
operators of positive chirality and one of negative chirality, which in this case has to be the third 
one due to the just mentioned condition on rij. Thus, the expression (j4.1|) must be multiplied by 
two factors from ()3.21|) and one from ()3.22|) . As a result, it reduces to 

1 1 "l+"2 o 
fib 1 bR £ 

C(ot ni , dn-21 < ^ni+n 2 ) = 77 ■ (4-2) 

7r J o 

The 6 = 1 limit of this three-point function is well defined and gives 

C c= i(a ni ,a n2 ,a ni+n2 ) = — fi~ « . (4.3) 

7T' 3 

One can show that for other choices of chiralities and for generic central charge the result differs 
from (J4.2)) only by a power of 6 and hence it leads to the same expression (J4.3)) in the c = 1 limit. 
This shows that our definition of the c = 1 limit does not depend on the artifact of the choice of 
chiralities. 

As explained above, the two-point function is obtained from the three-point correlator by putting 
one of the Liouville momenta equal to 6 (n = 0) and integrating with respect to (minus) the 
cosmological constant. For (J4.3|) this procedure gives 

1 R 

D c=1 (a ) = r log/2, D c=x (a n ) = r-/^, n>l. (4.4) 

7i n 

Of course, the same results can be found from the explicit expression ()2.10j) for the two-point 
function at generic central charge. 

Notice also that, if we had used the reflection amplitude D^ e ^{a) given in (|2.9|) as a starting 
point, instead of the two-point function derived from the three-point correlator, the computation 
in this subsection would have yielded a divergent result for _D c= i(a ). This provides a confirmation 
of the fact that the two-point function that we use is the one which is relevant for c = 1 string 
theory. 

To conclude, we see that the renormalization (J3.2U)) is sufficient to produce finite results for 
sphere correlators in the 6 = 1 limit. The result is analytic in the compactification radius R and 
there are no additional singularities arising at the self-dual point. 
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5 Disk correlation functions in the c = 



1 limit 



5.1 Divergent contributions as contact terms 

Let us now turn to the correlation functions characterizing the theory on the disk with FZZT 
boundary conditions. In our analysis we will encounter a new phenomenon. It turns out that 
all correlation functions in the boundary theory at the self-dual compactification radius R = 1, 
even after the multiplicative renormalization ()3.2()j) of the bulk and boundary vertex operators, are 
divergent in the c = 1 limit. However, we are now going to argue that the divergent contributions 
have the interpretation of contact terms, namely degenerate world sheet configurations, and should 
be omitted from the physical answers. 3 In order to understand the nature of the divergences, it is 
useful to anticipate here the general structure of the correlation functions in the c = 1 limit, which 
will be studied in full detail in the following subsections. 

First, it is very easy to characterize the correlation functions of operators with only vanishing 
boundary momenta of the matter part. As we will see, in this case the leading term diverges as 
1/e and is an entire (polynomial) function of the boundary cosmological constants. Such terms 
are well known in the matrix approach to non-critical strings where they are called non-universal 
contributions and usually neglected. This is justified by the fact that they vanish after a finite 
number of differentiations with respect to fi B , which means that such contributions disappear 
from the correlation functions with enough number of insertions of the boundary cosmological 
constant operators. In addition, since such terms are entire functions of \i B , they do not contribute 
to correlation functions with ZZ boundary conditions, which can be obtained from non-trivial 
monodromy properties of correlators with FZZT conditions [11,33,45]. 

The meaning of these non-universal contributions becomes clear when one passes to the length 
representation of correlation functions [47,48]. Suppose we have a correlation function C(/j, b ) on 
the disk with FZZT boundary conditions described by the boundary cosmological constant fi B . The 
parameter fi B is conjugated to the length of the boundary, and the correlation function can then 
be thought as a Laplace transform of the corresponding correlation function with a fixed boundary 
length 

oo 

C(ii b ) = J d£e~ e ^C{£). (5.1) 
o 

Let us concentrate on a contribution to C(fi B ) which is polynomial in fi B . Applying the inverse 
Laplace transform, one finds that it contributes to C(£) only at I = 0: 

ioo 

J dfi B e-^B p^ b ) „ p(-d/d£) 6{£) . (5.2) 

—ioo 

Thus, the contributions which are entire functions of the boundary cosmological constant corre- 
spond to contributions of degenerate world sheets with boundaries shrunken to a point, as shown 
in fig. n The same conclusion remains true if different segments of the boundary are characterized 
by different boundary cosmological constants. Such degenerate configurations should be excluded 
from the sum over two-dimensional geometries and, correspondingly, from the final result. 4 

3 The picture presented here arose in discussions with Ivan Rostov. 

4 These degenerate contributions may in fact still contain some physical information. Moreover, in some cases 
they do represent the physical part of the correlation functions as it happens, for example, in the case of topological 
gravity corresponding to c = — 2 [14]. However, the c = — 2 case is special due to the existence of two different 
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Figure 1: Degenerate world sheet configuration with shrinking boundary. It looks like a sphere with a 
puncture appearing at the place of the degenerating hole. 




Figure 2: The loop gas representation of correlation functions in the SOS model, a) A typical contribution 
to the sum over triangulations. b) The loop gas representation of a boundary operator (n = 2). In the 
continuum limit all lines come from the same single point of the boundary. 



The case when a correlation function involves boundary operators with non- vanishing boundary 
momenta is more complicated. In this case there are divergent terms up to order e - ™' '" 1 , where 
n-tot — ^2 n i is the total momentum of the boundary operators. The leading divergent terms are 
still polynomials in the boundary cosmological constants. It is these leading contributions that 
were calculated in [32], although our answers even for these terms are different from those of [32]. 
The reason for this discrepancy is that the c = 1 limit requires to be carefully denned and cannot 
be taken without considering the consistency of the coupling between the matter part and the 
Liouville part. 

Subleading divergent terms are not entire functions anymore. Therefore, the previous reasoning 
based on (|5.2j) . which allowed us to neglect the divergent contributions, cannot be directly applied 
and some additional argumentation is needed. A key observation is that the boundary correlation 
functions with momenta (j3.13|) can be written as polynomials in the one-point function Wis) of 
the boundary cosmological constant operator, which was given in (|2.25|) . This fact implies that 
all divergent terms come from the divergences of this function, and therefore one can remove all 
divergent contributions by renormalizing W(s). 

To understand the origin of this renormalization, it is useful to consider how these correlation 
functions are represented in the discrete approach to 2D quantum gravity. In this approach the 

theories with the same matter central charge: the above mentioned topological gravity and the so called bosonic 
string embedded in —2 dimensions [49]. Their correlation functions can be obtained either as finite parts of the 
limits b 2 — > 2 and b 2 — > 1/2 respectively, or as infinite and finite parts of the latter limit. However, we believe that 
two-dimensional or c = 1 string theory is defined unambiguously and its boundary correlators are given by finite 
non-degenerate contributions. 
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Figure 3: The two-point function d n (t\,l2)- 




Figure 4: Possible degenerations of the two-point function. The first picture shows the case when the 
boundary of an external piece shrinks, whereas the second represents a degeneration of an internal sector. 
In both cases, the degenerate piece becomes a sphere growing from the marked point. 

system we are dealing with, namely Liouville theory coupled with a matter CFT, can be realized as 
the so called solid-on-solid (SOS) model on a fluctuating lattice. 5 We are not going to present the 
SOS model here, and we refer the reader to [39] for the relevant details. What is important for us 
is the loop representation of this model [49,51,52], according to which the matter fields can be seen 
as non-intersecting loops on randomly triangulated surfaces. Correlation functions are given by a 
sum over all triangulations and all admissible configurations of the loops. A typical configuration 
contributing to the sum is shown in fig. |2K- 

The degenerate boundary operators with momentum p n = n/b in the matter part are realized 
in the loop representation by the so called "star" operators, that create 2n lines at some point of 
the boundary (see fig. |2t>). These lines are also non-intersecting and end at another point of the 
boundary. This fact allows to represent the general structure of the boundary correlation functions 
and was the starting point for their derivation from the discrete approach [39,53]. 

Let us then study in this context the structure of the correlation functions we are going to 
compute. The two-point function d n (Ti,T2) of boundary operators with momenta p n is presented 
in fig. El One can see that d n is composed of the product of 2n two-point functions do of operators 
with vanishing boundary momentum. The precise relation is easy to write in the ^-representation 

5 The corresponding matrix model realization was constructed in [50]. 
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Figure 6: Possible degenerations of the three-point function. The two cases correspond to shrinking of 
the boundary of an external or an internal sector, respectively. The degenerate piece appears as a sphere 
attached to the point at the boundary where two of the three boundary operators collide. 



due to the factorization property of the path integral measure [51-53] 

oo 

d n (£i, £2) = J < ■ ■ • d£' 2n do(£i,£[) d (£[, £' 2 ) ■ ■ ■ d (£' 2n , £ 2 ). (5.3) 

Q 

From this representation one sees that singular contributions to do give rise to singularities of d n , 
and the latter can have a much more complicated form. However, we know that the divergent 
contributions to do arise from degenerate world sheets. In fact, since d n can be represented as a 
polynomial in W(t), all its singularities can be interpreted in terms of degenerate configurations of 
some of the 2n + 1 disks composing the initial world sheet. The possible degenerate configurations 
are shown on fig. and correspond either to a disk with two coinciding marked points or to 
a pinched disk, together with a sphere growing from the marked point at the boundary. These 
contributions have again the interpretation of contact terms and should be neglected. Our analysis 
suggests that, in order to eliminate all of these contact terms, it is sufficient to remove the ones 
appearing in the correlation functions with vanishing momenta corresponding to each portion of 
the disk in fig. El We will see in subsection 15.41 that this is indeed the case. 

A similar analysis can be applied to the case of the three-point function C nijn2in3 (ri, r 2 , r 3 ). The 
representation in terms of non-intersecting lines is shown in fig. 0D We again see that the initial 
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Figure 7: The bulk-boundary two-point function. 



correlation function is reduced to a product of two- and three-point correlators with vanishing 
momenta. This is the origin of the representation of C ni)n2jn3 in terms of W{r) that we will find 
below, and allows to interpret each divergent contribution as a result of the degeneration of some 
of the disks which divide the world sheet in fig. The different possibilities are presented in fig. H 
and correspond to the same situations which we found in the case of the two-point function. Notice 
that one could expect the appearance of the degenerate configuration with all three marked points 
coinciding, which would give a world sheet in the form of a trefoil. However, one can show that the 
corresponding divergent term in Co,o,o vanishes. 

One can finally analyze the bulk-boundary correlator in the same way. Its loop representation 
is depicted in fig. [7| In this case the lines corresponding to the boundary operator, instead of 
connecting two points of the boundary, end at the point where the bulk operator is inserted. The 
naive counting would then tell that the correlation function should diverge as e~ n_m_1 , but we 
will see in subsection 15. HI that only the 1/e contribution is present. This can be explained from 
the fact, a consequence of (j5.31|) and (j5.12|) . that all terms divergent as e~ k with k > 1 in the 
boundary two-point function are proportional to (/! s (ti) — /i B (r 2 )) fe ~ 1 and, therefore, vanish at 
A*s( r i) = Vb( t 2)- Since the bulk operator cannot change the boundary conditions, the boundary 
cosmological constant indeed remains the same along the two boundaries of the disk formed by 
the outer lines coming from the two punctures. Thus, all contributions with a higher degree of 
divergence than 1/e vanish. The remaining divergent contribution turns out to be a polynomial in 
fi B , so it has the usual interpretation as contact term and has to be subtracted from the physical 
answer. 

After we remove all divergent contributions which can be interpreted as contact terms, all 
correlation functions become finite. This is a non-trivial fact which shows the consistency of the 
proposed renormalization. We interpret the resulting finite expressions as physical results for the 
correlation functions in Liouville theory at cl = 25. Notice that, in the cases where a comparison is 
possible, they also agree with the corresponding results of the discrete approach obtained in [39,53]. 

5.2 Bulk one-point function 

First, we compute the simplest correlation function on the disk, namely the one-point bulk corre- 
lator ()2.16|) . 6 Substituting the expression for the Liouville momentum from (j3.10j) and replacing s 

6 Some of the results of this subsection can be found in [45] where it is also shown that they agree with the results 
derived from matrix quantum mechanics. 
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by r using (|3.18|) . one obtains 

U(a n \r) = 1 [^ 7 (6 2 )]K^^)r (6 2 - f ) T (1 - £ - &) cosh ((l - £ - r) . (5.4) 

Taking the limit 6 — > 1 and using the renormalization (|3.21|) corresponding to the case of positive 
chirality, one finds, for n > 1 and generic compactification radius R, 



U c=1 (a n \r) = lim 

b-»l 



1 



A H^-^Jr (1 - 6 2 + f ) r (1 - j, - &) cosh ((1 - j, - r) 



7T 2 6 

^ n 



7rn nn 

R 6111 _R 



cosh Q . 



(5.5) 



The cases i? = 1 and n = should be treated separately. This can already be seen from the fact 
that the result ()5.5|) diverges at these values of the parameters. Let us start from the case of the 
self-dual radius with n ^ 0. The corresponding renormalized correlation function can be obtained 
from the first line in ()5.5|) substituting R = b. This gives 



U{a n \r) 



1 1 ( 71+1 



TP 



•Op ( n + 1 - b 2 ) r (1 - Bfi) cosh ((1 - sfi) r) 



(5.6) 



Substituting b 



e and expanding in e, one finds 



U(a n 



T 



27T 2 /T, 



1 



[n + l)e 



+ log/t + Cu(n) I cosh(nr) + 2rsinh(nr) 



+ 0(e), (5.7) 



where C[/(n) is some r-independent function of n which we are not interested in. We observe that, 
despite the renormalization, the correlation function still diverges. However, since cosh(nr) is a 
polynomial in coshr, the divergent term is a polynomial in the boundary cosmological constant 
p, B . As we explained in the previous subsection, such terms are interpreted as contact terms and 
correspond to degenerate world sheets. We then ignore them and take the "physical" one-point 
bulk correlation function in the c = 1 limit to be just the last term in ()5.7)1 7 



U. 



c=l 1 



\n+l 



a„.r 



ir 2 n 



■ fi 2 t sinh(nr). 



(5.* 



However, the case n = still remains singular. To get the correlation function in this situation, 
we return to the representation (|5.6jh put n = and expand around 6 = 1. The result reads 



U(a \ 



47T% 2 



Arc 2 e 



c 2 + c 3 log/2 



log 2 // 



+ 0(e) 



(5.9) 



7 Another possibility would be to find some "natural" subtraction defined before the limit b = 1 is taken. For 
example, one can define the renormalized correlation function as follows 

(i) (K 



U {Tcn) (a n \T) = U(a n \r) - T(l + b 2 - 2ba n )T(2a n /b - 1/b 2 - 1)P« ^ 



where P T l 1 ^(cos(a;)) = cos(nx) is the Chebyshev polynomial of the first kind. Then in the c = 1 limit one obtains 



U^ 1 \a n \r) = ]hnU^(a n \T) 



(-1)" + 1 A n 



fi 2 (log/2 cosh(nr) + 2rsinh(nr)) 

b— >l zn-'n 



This result is finite but still contains a non-universal piece. Therefore, we prefer to use the definition (|5.8|) . 
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Cj being some unimportant constants. The only term that depends on the boundary cosmological 
constant in a non-trivial way is the last one in the finite part of ()5.9|) . Therefore, in this case the 
"physical" correlation function is 

U c=1 (a \r) = (5.10) 

Z7T Z 

Note that it does not coincide with the naive n — > limit of ()5.8j) , but it agrees with the finite part 
of the n — > limit of (|5.5j) . 

Let us also give here the result for the one-point correlation function of the boundary cosmo- 
logical constant operator. In this case we have to use the multiplicative renormalization (|3.22|) 
corresponding to negative chirality 8 , and the correlator (J2.25j) becomes 

2b . 1 cosh(r/6 2 ) 

W{f) = p& ■ / J - 5.11 

7i sm(7r/tr) 

The substitution b = 1 — e leads to the following expansion 

W{t) » (~ - | +log/i J fi B + : ^v / A rsinhr + (5.12) 

so that we can take 



2 . £ g or , r B 2 



W c=1 {t) = ^-^/fLTsmhr. (5.13) 



n 2 



This is the well known expression for the one-point boundary function of c = 1 string theory which 
appears also as the resolvent in the matrix model formulation. Thus, we confirm the expectation 
that one should disregard the contact terms in order to reproduce the physical results of the c = 1 
theory. Of course, the universal contribution ()5.13j) agrees (up to contact terms) with ()5.1()j) . 

All correlation functions we are going to consider in the following are simplified only if the 
restriction R = 1 is imposed. Therefore, from now on we will concentrate on this particular choice 
of the compactification radius and omit the label U (R = 1)". 



5.3 Bulk-boundary two-point function 

Let us then proceed with the analysis of the bulk-boundary correlation function. Substituting 
(J3.13j) into (J2.18|) . one finds the following expression for the bulk-boundary correlator 

R(a n ,p m \r)= ^ m {b 2 )b 2 - 2h2 Y^ KtKz, (5.14a) 
1 n(b+l)T b (^-b)T b (b-f)T b ("-±l-b)T b (2b-r t y ^ b) 



J S 6 (t + ^f±^)S b (t + 6+^±i) 



(5.14c) 



8r 



The choice of chirality can be explained as follows. W is obtained from the boundary three-point function by 
integrating twice with respect to fi B . The three-point function, according to our convention, has two operators with 
positive and one with negative chirality. Since each integral corresponds to removing one boundary cosmological 
constant operator, which is an operator of positive chirality, the resulting one-point function then necessarily corre- 
sponds to an operator of negative chirality. This implies that, strictly speaking, W is not the one-point function of 
the boundary cosmological constant operator, but rather of the operator B- q = e tqx e b ^ . However, this reduces to 
Bq in the c = 1 limit, so we keep referring to W using this simpler terminology. 

Another related subtlety is that in the theory with background charge the disk partition function is not well 
defined since the background momentum remains non-compensated. 
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The second factor IZi can be simplified by using the recursion relations QA.5[) . The result is 

2m 

(2n) 2m b n+m+1 - 2b2 - m l^ ±1 m\T(b 2 -n) U F (*) 

n x = — — — . (5.15) 

mm m~n—l m+n+1 *> ' 

r(2m + i-^ )n r(A) n r(i + A) n r(i + £) n r(£-i) 

k=l k=—m k=~n k=n+l 

The main non-trivial problem is to compute the integral 7^2- Thanks to the recursion relations 
among S& functions, the integrand can be simplified to just a product of sines: 

+ioo 

K2 = ^Jr~ I dt^^ . (5.16) 

ft Ml (*" (*+{))] 

I n— m+1 

K ~ 2 



22m+l 

— ioo 



The definition of the contour of integration in the t-plane is the one explained in section 12.21 the 
contour lies on the right of the poles of the integrand coming from the numerator, and to the left of 
the poles coming from the denominator. We also mentioned that this definition of the contour can 
become problematic if poles of the numerator become coincident with zeroes of the denominator. 
However, this never occurs if we limit ourselves to the case n > m, so in the following we assume 
this restriction. Closing the contour in the right half plane, the integral is then computed by picking 
residues at the following poles 

t = 3 - + lb, l>0 and t=- J T +lb, />0, (5.17) 

b b 

where j = n ~T ,+1 , . . . , m+ ™ + . The computation of the residues gives: 

-2rf A+0 V- „-27-f-&+0 



(5.18) 



2 2m . irU II sin (J - k)) n sin (f (j + k)) 

2 fc^j fe 

(-l) m tt6 >A - e^" 2T 

2^ l-e-V £^ nsin^-^nsin^Q' + fc))' 

2 k^j k 

where to avoid clutter of indices we have understood the range of the products over k, which is the 
same as in (|5.16|) . Putting (|5.15|) and (|5.18|) together, and taking into account the renormalization 
of the vertex operators (|3.20j) . we finally get 

. , . 2m 
(_l)m m! ^2m-l 6 -l^^ 2 ±l-l r(1 _ 6 2 + n)r (2^_ 1 ) jj r ( fc ) 

R(cn n , d m \r) = : fc-m+i 

v ' ' ' I mm to— n— 1 m+n+1 

r(2m + 1 - *») n r (£) n r(i + £) n r(i + £) n r(£-i) 

k=l k=—m k=— n k=n+l 

n+m+l „. „. (O.lyj 

1 2 — 32£ ^Z-2t 

I \ --, e b2 — e^- 



X 



l - e- 2T 



n+_m+l 



,_ n-m + l • 2 

■> 2 



sin^ II [ sin {§ (J - k)) sin (j + k))] 



n — m + 1 
" 2 
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Let us now extract the leading terms in the b — > 1 limit. We find 

R(a n ,(3 m \r) 



(5.20) 



\n+l 



n+m— 1 / 777, 



2vr% 



n 

K k=l 



T (k + m)T (k + n — m + 1) 
T(k)T(k + n) 



(l + c R (n, m)s)TZ njm (r) + 0(e), 



where 



, N ™ , sinh((l- 2 -1 + " +1 )t) 



3=0 



3=0 



sinh r 
sinh ((2j + n — m)r) 



sinh r 



+ 2(2j + n - m + l)e 



r cosh ((2j + n — m) r) 
sinh r 



(5.21) 



+ 0(5) 



with vAn, m) 



(-rp (j+n-m)! 



U]Vl , ,,<>) — .y ( m _j); (j +n+1 )i ■ F° r simplicity we also included the term proportional to log fx into 
the constant cr(ji, m). As a result, the leading term of (|5.2U|) is proportional to 



■^v i (n,m)P : 



(2) 



2j+n-m-l I ^/jjj 



(5.22) 



3=0 



where P^} 1 (cos(x)) = s ^?^ is the Chebyshev polynomial of the second kind. Thus, the leading 
term is divergent but is a polynomial in fi B : it can then be interpreted as a contact term and 
disregarded. The physical contribution can be read off the second term in (|5.21|) and is given by 



R c=1 (a n , /3 m \r) = — fi 



m\ „ n+m-i / t— r r (k + m) r (k + n — m + 1) 



7T^ 

m 

j=0 



n 

^=1 



r(fc)T(A; + n) 



(— !)"+■? (2j + n — m + l)(j + n — m)! r cosh ((2j + n — m)r) 



(5.23) 



(m — j)! (n + j + 



sinhr 



In the case where the condition on the momenta following from the matter sector is taken into 
account, requiring n = m, this expression gets some additional simplifications, becoming 



Rc=i(a n , p n \r) 



. I n 



7T 



E 

3=0 



_!)n+i (2j + l) ( n !)2 T cosh (2j T ) 



(n — j)\ (n + J + 1)! sinhr 



(5.24) 



Let us also perform a consistency check, namely that the result ()5.20j) together with (|5.21|) 
reproduces the bulk one-point function found in (|5.7|) for n > and in (J5.9|) for n = 0. For m = 
one finds 



R(a n ,/3 Q \r) 



2tt 2 



sinh (nr) , , 2r cosh (nr) 

1 + c R (ra, 0)e) + 



[n + l)e sinh r 



sinh r 



(5.25) 



Then the integral with respect to —fi B is easily calculated, yielding 



J dfi B R( 



a n ,/3 \r) 



2ir 2 n 
2n 2 



1 



(n+l)e 



Cr(ti) ) cosh (nr) + 2rsinh (nr) 



n > 0, 
n = 0, 



(5.26) 



in perfect agreement with ()5.7j) and with the universal part of (|5.9j) . The non-universal part for 
n = does not depend on r and cannot be fixed starting from the integral (|5.2fij) . 
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5.4 Boundary two-point function 

The next function which we are going to consider is the boundary two-point correlator. Although it 
follows from the three-point function, it is useful to start with this simpler quantity since it allows 
to illustrate the main non-trivial points of the renormalization leading to a finite answer in the 
limit b — *■ 1. 

The explicit expression for the two-point function was found in appendix and is given in 
(|2.24jl together with ()2.23|) . Using (J3.13J) and replacing by 7$, one obtains 

1 (2n+l_ 1 \ 

2(^±l-b) \n^(b*) b^ 2 } ^ J T b (b- T- 1 - b) 

d(f3 n \ri,T 2 ) 



s 6 (&-! + ^) s 6 (&-f-i^p) s 6 (&-f + i^p) s b (&-f-i^p)- 

(5.27) 

By using the properties of the special functions, the factor in the numerator can be transformed as 
follows 

q-g^-O Mfc-y) 2 ff 1 

r 6 (^ti _ r (2n + 1 - b 2 ) r (^±±) Al 2 

For the denominator in ()5.27j) . we get 



sin ^ ' 



(5.2* 



n 

2 4n FT Tsinh ( n+Vi+2mk \ ; U ( T 1 -T 2 +27nfc \] 



S 6 (6 - ? + S 6 (& - f + i^f) sinh (a±a) sin h (a=a) 



(5.29) 



Then the final expression, without special functions and renormalized accordingly to ([3.20)1 . turns 
out to be 

2n+l n 

rf(fi,KT 2 ) = ^K^-') *- g-t=^ . (5.30) 

7T cosh ri — cosh r 2 

Remarkably, this expression can be rewritten in terms of the (renormalized) one-point functions 
Wi = Wiji), i = 1, 2, given in (j5.11|) . The result reads 

, s W,-W 2 tt / - 9 -9 2fc7T ^ ^ 46 2 i sin 2 ^\ 

diPnln, r 2 ) = -a n -A 2 - T] K + Wf - 2 cos — - WiW 2 - — -jf- , 5.31 

fi 12 y " 6 2 tH sm y 

where we used the definition ()3.19j) and introduced 

( 7T \ 2n -r-r Sin £ , 

In particular, for n = 0, namely for the two-point function of £>g lT2 and B 1 }^ 1 , we obtain 9 

_ 

d(b\n,T 2 ) = — 1. (5.33) 

Ml2 



9 In the limit t\ = T2, (|5.33|) reproduces the result ljC.19fl for the correlator of two cosmological constant operators 
multiplied by the renormalization factor -tT(1 — b 2 )T (l/b 2 — l). 
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Note that the expressions (|5.31|) and (|5.33|) . with a slightly different normalization, appeared in [39] 
where they were derived from the discrete approach based on the SOS model and the loop gas 
representation. 

To extract the physical result for the c = 1 theory, let us note that all coefficients in (|5.31|) in 
this limit remain finite. In particular 

lima n = ^— - rr . (5.34) 

Thus, to get a finite result, it is enough to renormalize the one-point functions W{t) appearing 
in ()5.31|) . The interpretation of such a renormalization was given in section 15.11 where it was 
shown that it corresponds to neglecting world sheet configurations with two marked points on the 
boundary close to each other so that the world sheet degenerates. Replacing then W(r) by its 
physical part in the c = 1 limit given in (J5.13j) , one finds 

. , ai x (-lr+V* 1 W c=l (n) - W^ x {t 2 ) -A- / 2 iQk 2 fi\ 

d^n\n,r 2 ) = 22n(2n + 1)] Ab(ti) _ Ab(t2) n (0»^fr) - W ~*M) --^) 

= ' ( :^ r+1/i ^ riSin !! ri " r2S : nhr2 ff Undnhn - r 2 sinhr 2 ) 2 - (2nk) 2 ) . (5.35) 
7T 2n + 2 (2n + 1)! cosh n - cosh r 2 11 U ^ 



It is clear that the divergent contribution to W(r) gives rise to a series of divergent contributions 
to d((3 n \r 1 ,r 2 ). Only the leading contribution is an entire function of fi B , whereas all others have 
a non-trivial functional dependence on Tj. Nevertheless, all of them correspond to contact terms 
according to our reasoning in section 15.11 The only physical contribution is given by the finite 
expression (|5.35|) . 



5.5 Boundary three-point function 

Now we turn to the most complicated correlation function, the three-point correlator ()2.20|) . Sub- 
stituting the momenta ()3.13|) and replacing Sj by n, one obtains 



Ci 

Co 



l2\ L2-26 2 



i+E n 

2b' 1 



C\ C 2 C3 



Tb + »i-»2-n3 ) Tb (_ b + "i+" 2 +"3+2 ) ^ ( "i+" 2 -"3+l ) ^ ^ 



"2+"3 + l N 



r b (6 + i) r b (-6 + ^±i)r 6 (-6 



s.(^ + ^) s 6 ( 



W3+1 



— I 



; ri+T3 



27rft 



) 



A T 2 +T 3 \ 

1 2nb ) 



-Aici 



dt ] [ 



S fc (t + Uj) 
~S 6 (t + V;)' 



(5.36a) 
(5.36b) 

(5.36c) 
(5.36d) 



and 



u 3 
u 4 



"1+1 _|_ a 21+22. 

ft TJ 2!ri) ' 

"1+1 _ A T1-T2 

ft 27Tft ' 

b-^ + i T -^ 



"2 + 1 
ft 



2?rft ' 



_|_ A T1-TZ 



"l+"3+2 
6 

"1~"3 + 1 

6 

1 



* 2tt6 ' 



2-Kb ' 



b+l 



(5.37) 



Let us start by simplifying C\. In the following we will restrict ourselves to the case when 
the momenta rti satisfy the triangle inequality, \n 2 — ns\ < n\ < n 2 + n 3 . This allows to avoid 
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considering many different cases which can not be treated simultaneously and is consistent with 
the momentum conservation imposed by the matter sector. One then obtains the following result 
for the quantity in (|5.36b|) : 



n2+n3— ni 2n 2 2n 3 



,3+3n 2 +3n 3 -n 1 -2b 2 -£^±i II T i 1 ~ p) II T (£) II T (£) 

^7 ^=0 fc=ni+n2-ri3+l fc=rii+ri3-ri2+l gg\ 

y ] r(2n 2 + i-6 2 )r(2n 3 + i-6 2 ) n r (-1 + £) 

fc=2ni+l 

Now we pass to the evaluation of C3. As a first step, we can eliminate four of the eight 
functions appearing in the integral, obtaining 



Til 



C, =— ^ / (it 



— too 

Vti+n.3+1 ni— n.3 

n «in(f (t + ^ + f)) n sin(f (t + ^ + f)) 

. fc=n 2 +l k=— 712 



(5.39) 



The contour of integration in the complex t-plane is defined as usual to lie on the right (left) of the 
poles of the integrand coming from poles (zeros) of the S& functions in the numerator (denominator) 
of (j5.36d|) . Due to our condition on the momenta, the poles coming from the numerator and from 
the denominator never come to collide, and the contour is well defined. 

However, there is an important difference between this integral and the one appearing in the 
expression for the bulk-boundary correlator. Whereas in the latter case we could close the integra- 
tion contour in the right half plane due to the exponential fall off of the integrand, this cannot be 
done in the present case. In a similar situation, the one of the derivation of the boundary two-point 
function from the three-point function in appendix as a solution to this problem we proposed to 
use some "shift equations" characterizing the result of integration. We apply the same strategy in 
appendix ID. II for the case at hand. We demonstrate that the equations we derive for the integral 
()5.39|) are satisfied by the expression that one obtains by closing the contour of integration in the 
right half plane and taking into account only contributions from the zeros of the sine functions in 
(|5.39j) . Notice that, according to these "effective integration rules", the poles of the S^-functions, 
which have a much more complicated form, remarkably do not contribute to the result. 

In summary, in order to compute the integral we have to pick residues at the following two 
series of single poles 

i T^n_J_ b fi = n a + l,... ) n 1 + na + l > I > 0, ( } 

3 2nb b ' \ J = -n 2 ,...,m-n 3 , I > 0. K J 

For both series the residues of the integrand 1{t) can be written as 

Res t=tj , Z(t) = — 1 , J x S jh (5.41) 

4"i J! sm \w (k- 3)) 

where k in the product runs through the union of two intervals k G [— n 2 , ni—n 3 ]U[n 2 + l, 71.1 + 723 + 1] 
and the quantity Sji is defined as 

31 s 6 ^ + v + (^ + i)&)s fc H I l^ + ¥ + ^ + 1 ) & )' 
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All of the special functions in Sji can be eliminated by multiplying (|5.42j) by (|5.36c|) . getting 



(5.43) 



where 



ni-j 



Aj= EI [sinh ( T1+T |^ 2 " fc ) sinh ( 



k=-n 3 



Ti — T°3-\-2nik 
2b 2 



112 



)] n [sinh ( a±^3gg ) sinh { T2 - T £ 2nik )} , (5.44) 
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l-l 

Y\ [sinh ( 

fc=0 



n+T3 



Til 



k=l-j 

kb 2 ) sinhf^^p + vr^ 2 )] 



sinh (a±s) sinh (^p) f] [sinh ( 



T2+T3 



7T2 



k=0 



kb 2 ) sinh + nikb 2 )] 



(5.45) 



Collecting all factors from (j5.41|) - (j5.45J) . we can write the final answer for the product C 2 C 3 : 

00 ni+713+1 



C 2 C 3 = 27ib{-l)^ n H n2+ns - 1 < J2 B i m 



(✓4-7 



+"3 + 1 



;=1 j=n 2 +l 



00 ni—n-j 

+E* E » 



fc— n<2 + 1 



ni-n 3 

n sin n sin^^-j)) 



(5.46) 



+n 3 +l 



ni-n 3 



<=o i=- 2 [I sin(^(fc-j)) EI sin(£(fc-j)) 

A) = ?l2 + 1 fc = -n 2 

It turns out that the sums over / in ()5.46|) can be calculated explicitly. For this let us introduce 
the shift operator in r 3 

A 3 = e- 2mb2a ^. (5.47) 

Using this operator, one can write 
4^ 



1=0 



El ^ cosh n— cosh T3 

V ^ cosh T2 — cosh (T3 — 2-7T «6 2 ) 



4^ 



^ CQsh Tl —COsh T3 

^ coshr2— cosh(r3 — 2irib 2 ) 



1=1 



fi2z cosh r 2 — cosh (r 3 — 2nib 2 ) /2 23 coshr 2 — cosh(r 3 — 2nib 2 ) 



(5.45 



Then from the following identity 

cosh Ti — cosh T3 



1 - A, 



it follows that 



coshr 2 — cosh(T3 — 2nib 2 ) 
4/2 



coshr 2 — cosh(r 3 — 2ixib 2 
cosh T\ — cosh r 2 

4/2 



-1 



^12 ^23 



E« 



^12 ^13 



(5.49) 



(5.50) 



1=1 r ^ r ^ l=Q 

Remarkably, the results do not depend on b. This means that they will not generate an expansion 
in e = 1 — b. 

Note that the quantities Aj defined in ()5.44j) can be written in terms of the renormalized one- 
point functions W(r) given in (|5.11|) : 



Abfil 



7r sin(7r/6 2 ) 



I] (W(n)-W(7s-27riib)) EI (W(r 2 ) - W(r 3 - 27riA;) 

fc=— 713 ^ ' k=l—j ^ 



(5.51) 



25 



Thus, after the multiplicative renormalization (|3.2(J|) . where we choose the third operator to be the 
one with negative chirality, the expression for the full three-point function reads as follows 

C ( P n ! , P n2 , P n 3 | T\ , T 2 , T 3 ) 

( , n 1+ n 3+ l II (W(n) - W(r 8 - 2^)) EI (W(t 2 ) - W(t 3 - 2*1*) 
I 1 v fc=— rt3 ^ ' k=l—j ^ 

° I Al2 ^23 ^ ni+n 3 +l . m-na 



+ 



ttt; E 



I! sta(£(*-j)) n sin(£(fc-j)) 

fe=n2+l k= — U2 

, II (W(n) - W(t 3 - 2™*)) n _ (W(r a ) - W(t 8 - 2*1*)) 



k=—ri3 v 7 fe=l— j 

A*12 A*13 



J=-n 2 



711+713+1 7ll— "3 

EI sin(£(*-j)) II sin(£(fc-j)) 

fc = ?l 2 + l fc=-n 2 



where the overall coefficient is 



r(2 Wl + 1 - 5 2 ) T(2n 2 + 1 - 5 2 ) r - 1) / n x E 
C ° l Sm ¥) Cl (5 ' 53) 

"■2+713— "1 2n 2 2n 3 

r(2m + 1 - & 2 )r (a$ti - 1) n r(i-£) n r(£) n r(£) 

fc=l fc=ni+n2— 713+1 fc=ni+n 3 — n 2 + l 



_v^„_n ni+n 2 +n3+l 
2 En 17 r-2ni63n 1 -n 2 -n3 ( s i n iL) ^ 1 r(2n 3 + 1 - 6 2 ) f] F (~1 + 

fe=2m+l 



Although this result is still lacking of the explicit cyclic symmetry under simultaneous permutations 
of rii and Tj, we checked on examples for low values of momenta that the cyclic symmetry is indeed 
present. Several explicit expressions can be found in appendix ID. 21 From these examples one also 
observes that all imaginary terms produced by shifts of the argument of W in (|5.52|) are canceled 
and the final expressions can be written solely in terms of Wfa). Actually, this is a consequence 
of the cyclic symmetry, since the imaginary terms are odd functions of r 3 , whereas the function 
(J5.52j) is even in T\ and r 2 , which means that it should also be even in r 3 . This is a crucial result 
for the c = 1 limit. Unfortunately, we cannot prove it in general explicitly due to the complicated 
combinatorics involved in (|5.52|) . but since it follows from the fundamental symmetry of the three- 
point function, it ought to be true. 

The general structure of the three-point function evaluated for degenerate momenta, namely 
being a polynomial in W{ji) divided by the differences of the boundary cosmological constants, was 
derived from the microscopic approach in [39]. However, the explicit expression for the polynomial 
was missing. Although our result is also not as explicit as one might wish, in order to obtain this 
polynomial it is enough to expand W{jz — 2irik), cancel all imaginary terms and rewrite the result 
in terms of W{t 3 ), as we did for some cases in appendix ID.2I 

The fact that the three-point function is a polynomial in W{ti) allows to apply the renormal- 
ization procedure described in section I5~T1 and already successfully used for the boundary two-point 
function. Let us first study the limit b — > 1 of various coefficients appearing in (|5.52j) . Expanding 
b = 1 — e, the leading part of the factor Cq given in (|5.53|) reads 

2n 2 -l 2n 3 -l 

n k\ n v. 

C ~ (27TC) 2W1+1 ( l)|("2+7i3-ni)(n 2 +n 3 -ni~l)+l /7T\^ n ' fc=rai+n 2 -n 3 fc=ni+n 3 -n 2 (5 54) 

~ ^ ' ^ ' \l) rii+ri2+n 3 —l n2+n 3 —ni ' \ ' / 

II k\ U k\ 

k=2ni k=l 
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Next, the coefficients of the polynomials in W{ti) appearing in (|5.52|) become 

1 _ (_l)«3+i+l ( n;j _ ni _ I + J )| ( 2vre )-2n 1 -l 

ni rf +1 • 7* a. ^ "tt 3 • (■« a ^ ~ K + j)' 5 ~n 2 - 1)! fai + ra 3 + 1 - j)! 
I] sm(^(A;-j)) n sm^(fc-j)) 

fc = n 2 + l k = — ri2 



(5.55a) 



1 ^ (_i)»i-t«.+i+i {n2 _ j}! ( 27r£ )-^-i 



ni tt +1 • ^"tt 8 • f*^ ^ («2 + i)!(ni-n 3 -i)!(ni+n 3 + l-i)! 

I] sm(^(fc-j)) n sm(^(A;-j)) 

Thus, the ^-factors from (|5.54jl and ()5.55j) cancel and the three-point function can be written in 
the following form 
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According to our previous observation about the cancellation of imaginary terms, the expression in 
curly brackets is a polynomial in W(ri). The coefficients of this polynomial are finite in the limit 
b — > 1. Therefore, all divergences come only from the one-point functions and can be interpreted 
as contact terms. To remove them, it is enough to replace all one-point functions by their physical 
expressions in the c = 1 limit, W / C=1 (rj), given in (J5.13|) . After this procedure, the resulting 
correlator, which we denote C c= i(/3 ni , (3 m , f3 n3 \ri, t 2 , t 3 ), will be finite and possess all necessary 
properties of a three-point correlation function. In particular, one can check that it reproduces the 
boundary two-point function found in the previous subsection. The details of this calculation are 
presented in appendix ID. 31 

6 Conclusions 

In this paper we have studied the correlation functions of Liouville theory in the limit in which 
the Liouville parameter b goes to 1 and the central charge cl goes to 25. The main motivation for 
considering this limit is that it corresponds to Liouville theory coupled to c = 1 matter, namely to 
two-dimensional string theory. 

Our primary interest has been in the particular case of two-dimensional string theory com- 
pactified on a circle at the self-dual radius R = 1. This case turns out to be interesting for two 
reasons. First, it allows a very explicit representation of the correlation functions, where all special 
functions characterizing the general results of Liouville theory disappear. Second, when the Liou- 
ville momenta are restricted to the values imposed by the compactification, additional divergences 
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appear besides the ones that are already present in the limit b — ► 1. Whereas the latter divergences 
can be removed by a multiplicative renormalization of the Liouville couplings, the former are more 
subtle, and some additional physical considerations were needed in order to define the theory at 
these values of parameters. 

Inspired by the picture emerging from the loop gas representation of non-critical strings, we 
found a physical interpretation for the divergences appearing at the self-dual radius. We associated 
them with various degenerate world sheet geometries, the simplest example of which is provided 
by a disk with vanishing boundary length. Such geometries must be excluded from the path 
integral and hence the corresponding terms, which we call contact terms, must be subtracted 
from the correlation functions. As a result, we proposed a precise procedure to take the c = 1 
limit based on the coupling of Liouville with a matter CFT representing a free scalar field with a 
background charge. The procedure goes as follows. After a multiplicative renormalization of the 
vertex operators, one evaluates the correlation functions for special values of momenta, takes the 
limit b — *■ 1, and finally subtracts all contact terms. 

In this way we arrived at a well defined set of finite correlation functions. The main results are 
the following: 

• three-point function on a sphere — (14.3(1 (for a generic compactification radius R) ; 

• bulk one-point function on a disk — (|5.5(1 (for generic R) and (|5.8(1 (for R — 1); the case 
of the cosmological constant operator is special and is given in (|5.1U() (for the bulk) and in 
(|5.13j) (for the boundary); 

• bulk-boundary two-point function — (|5.23j) (for R = 1); 

• boundary two-point function — (15.35(1 (for R = 1); 

• boundary three-point function — (15.56(1 (for R = 1), where in order to get the final result 
one has to perform an additional renormalization, by evaluating the polynomial in Wfa) and 
replacing the one-point functions by their c = 1 limit, W / C=1 (rj), given in (|5.13|) . 

In this paper we also found several new results relevant for Liouville theory at generic central 
charge. First, we derived the boundary two-point function from its definition in terms of the three- 
point function, which is the only sensible definition in the context of quantum gravity and string 
theory. As expected, it turned out to be proportional to the boundary reflection amplitude and 
the precise relation between them is given in (|2.24(l . 

Second, we obtained expressions for the correlation functions evaluated at the values of Liouville 
momenta corresponding to degenerate momenta in the matter part. For such momenta all special 
functions and integrals can be eliminated and the results are presented in a simple form. The 
expressions for the bulk-boundary, two-point and three-point boundary correlation functions can 
be found, respectively, in (|5.19(1 . (|5.31(1 and (|5.52(l . 

All of these results can in principle also be found by using the techniques of matrix models. The 
matrix model relevant for the c = 1 case is matrix quantum mechanics (MQM), which provides a 
solution to two-dimensional string theory both in the trivial and perturbed backgrounds [54, 55] 
(see [56] for a review). The comparison between MQM and CFT results for the one-point bulk 
correlation function on a disk was already performed in various papers devoted to the study of non- 
perturbative effects in two-dimensional strings. Essentially, most of the checks of the agreement 
between MQM and CFT explore one or another manifestation of this correlator. 

However, the correlation functions with non-trivial boundary operators are hard to consider 
due to a difficulty in introducing such operators into the MQM framework. Instead, they were 
investigated using statistical models on a fluctuating lattice and their loop gas representation 
[39,43,52,53]. The main result of this approach is the coincidence of some discrete equations for 
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the correlation functions (similar to the shift equations used in this work) with the corresponding 
equations satisfied by the CFT answers [39]. However, also these statistical models are typically 
formulated for generic values of the central charge, while the problem of the singularity of the c = 1 
limit remains. It would therefore be quite interesting to find a model defined directly for c = 1 
that is able to incorporate non-trivial boundary operators. 

Finally, let us remark that in the case of Liouville theory on a disk we only provided expressions 
for the correlation functions with FZZT boundary conditions. The corresponding correlators for 
ZZ branes [26] can be obtained by using a known relation between FZZT and ZZ boundary states. 
However, whereas a ZZ correlator in the bulk is expressed just as a difference of FZZT correlators 
evaluated at special values of boundary parameters [5] , the case when non-trivial boundary opera- 
tors are involved is more complicated [27, 29] . It would be nice to find a geometric interpretation 
for the corresponding relation, similar to the one that holds for the bulk correlators [33,45]. Our 
correlation functions might be useful for this purpose since they are much more explicit than the 
expressions at generic central charge. 
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A Special functions 

The special functions appearing in the main text are defined as follows 



T b (x) = exp 



dt 

T 



Sh(x) 



r b (x) = 



and we also use 



o 

r b {x) 



e -xt _ g-Qt/2 

;i - e- bt ) (1 - e-'/ b ) 



(x - Q/2f _ t x-Q/2 



T b (Q - x) 
1 



exp 



dt 

T 



smh((Q-2x)t) ^x-Q/2 



2sinh(6t) sinh(t/6) 



dt 

r b (x)r b (Q -x) exp \ J J 

o 



--x) e 



sinh 2 - x) | 
sinh |r sinh ^ 



¥ b n — 



dcf dT b (x) 



dx 



r b (b). 



x=0 



They possess the following properties 
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S b (Q/2) = 1, T b (Q/2) = 1, 

S 6 (x) = S^Q-s), T 6 (a;) = T b (Q-x), 

S b (x + b) = 2sm(nbx)S b (x), ? b (x + b) = 1 (bx)b 1 - 2bx T b (x), 

S b (x + l/b) = 2sm(nx/b)S b (x), T b (x + 1/6) = 7 (x/6)6 2 ^- 1 T fe (x) 



(A.6) 



r&(a;) is meromorphic with poles: x = —nb — m/b, n,m G Z-°; 
S b (x) is meromorphic with poles: x = —nb — m/b, n,m G Z-°, 

and zeros: x = Q + nb + m/b, n,m G Z-°; (A. 7) 

Tft(x) is entire analytic with zeros: x = —nb — m/b, n,m G Z-°, 

a? = Q + nb + m/6, n, m G Z-°. 

B Derivation of (7 from i? 

In this appendix we show how the one-point bulk correlation function on a disk can be derived 
from the bulk-boundary correlator. To perform this computation, one should first take the bound- 
ary momentum to represent the boundary cosmological constant operator, and then integrate the 
resulting expression with respect to —\i B - After the first step the expression (|2.18|) becomes 

Q-2a 1 

R{a,b\s)= 7c f i 1 {b 2 )b 2 - 2b2 2b ' (B.la) 

n = El (i)r 6 (2a-6)r 6 (j + 6-2a) 
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7 S t («-a + |+f.)S t (« + a+ A) 

— ioo 

The first factor can be simplified to 

1_ r(i-b 2 ) r(i + i-f)' 

whereas the second factor reads 



lOO 

g— 2nts 



1Z 2 = 7TI I dt 7 7 i^r-v 7 7 rrr. (B.3) 

J 2sin(f (t + «-i))sin(f (t-a + i)) 

— ICO 

According to the integration rules described in section 12. 2\ the integral is computed by picking 
residues at the following poles 

t = -a + ^- + lb, l>0 and t = a - ^- + lb , I > 0. (B.4) 

2b 2b 



The result is 
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oo oo 



2 sin (f (2a- I)) ^ 
nb sinh (7rs {2a — Q)) 



e —v" 2b 



sin (f (2a - J)) sinh (nsb) 



(B.5) 
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Taking (iFTTafl . (TET2j) and (lE~5j) together, one obtains 

Finally, the integral with respect to the boundary cosmological constant gives 

— / d/j,„ R(a,b\s) = —nb, — j- / ds sinh(7rs6) R(a, bis) (B.7) 
J V sm(7r6 2 ) J 

= \ [n^(b 2 )} {Q ~ 2a)l2b T(2ba - b 2 )T(2a/b - 1/b 2 - 1) cosh ((2a - Q)irs) . 
nb 

This reproduces precisely the correlation function (|2.16jh which differs by a factor 2n from the 
standard expression used in the literature. Its appearance can be ascribed to the residual rotation 
symmetry of the one-point function. 



C Derivation of d from C 

In this appendix we derive the boundary two-point function starting from the boundary three- 
point function. This is the only correct way to obtain it in the context of quantum gravity. The 
approach is the same as in all previous cases: to take one of the boundary operators to be the 
boundary cosmological constant operator and then integrate with respect to —\i B . However, since 
the expression (j2.20j) is not explicitly symmetric in three momenta, there are three possible ways 
to perform the calculation, all of which should give the same result. We are going to analyze all 
of the three possibilities, and this will give a non-trivial check of the symmetry of the boundary 
three-point function under cyclic permutations. 



C.l The case £| Sl = B s h s 

The first possibility is to take (3\ = b, j3 2 = P3 = P, and s\ = s 2 = s, S3 = s'. Then the three-point 
function can be represented as 



2\t2-2b 2 



" 3^ b r (i) n (i) r (i) 



u l u 2 L '. 
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C{b,p,p\s,s,s')= 7r/i 7 (6 2 )6 : 
r 6 (2/3-6)r 6 (| + 6-2/3)rg(|) Y(2pb-b 2 ) 
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r b (l + b)r b (l-b)r 2 (l + b-2p) ' r(i-6>) rfa + i-f) 
s b (l + b-p + i s -=f) s b (l + b-p- ijY) 
s b (p + i s -^-)s b (p-i s -¥) 

ioo 

S b (l+is+t)S b (l+t)s b (f3+i^+t)s b U+b-/3+i^+t) 



S b (2P), 



-Ani I dt 



Til 



S b (Z+b-P+i^+t)s b (l+(3+i^+t)s b (l+b+is+t)s b (l+b+t) 

dt 



4 J sm(7vbt) sink's + t)] sin [f (p + i^ + t)] sin [f (± - p + i^ + t)] 



(C.la) 
(C.lb) 

(C.lc) 



(C.ld) 



Note that, in contrast to the integral appearing in the expression for the bulk-boundary correlator, 
the integrand in (|C.ldJ) does not vanish when t — > ±00. Therefore, in principle, we cannot close 
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the integration contour either in the left or in the right half plane. Nevertheless, let us assume 
that it is possible and the contour is closed in the right half plane. We will justify this assumption 
shortly. Then, using the integration rules described in section the last factor can be evaluated 
by residues. This gives the following result 



fc=0 



sin(™6) I sin [ f (| +/3+ ^)] sin [ f (Mi_ /3+ ^)] sin [ f (| +/3 _^)] sin [ f (Ml_^_^)] 



+ 



sin[f (2/3— I)] I sin [irb(kb+l3-i^ )] sin \Kb{kb+l3+i^f )] sin[7rb((fe+l)6-/3-i^)] sin[7rb((fe+l)6-/3+i^)]' 



It is easy to convince oneself that this expression can be rewritten in the following form 



(C.2) 
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(1) = E 

2isinh(7rsfe) sin[f (2,(3-^)] 

X 



1-e "5" 



+ 



1-e 



^(cot[f(i-/3 +i ^)]-cot[f(^ + i^)]-cot[f(^ +i ^)] + cot[f(i^ + i^)]) (C.3) 
cot^rrb(b-(3+i^-^ -cotj^/M-i^!^)] -cot [wb^+i +cot Urbfb-fi+i s ^)] ) 



where we used a shift operator in s similar to the one in ()5.47|) . 
Let us consider the function defined by 



@h(x) 



1-e 



b / , \ 1/& / 

— rK- COtiTTOX) H j— COt(7TX/0). 

OOx 1 ~tO x 



1 — eb 



(C.4) 



Applying the operators 1 — e bdx and 1 — e* 8 * to this definition, one can establish the following 
properties under shifts of the argument 



@b(x + b) = @b(x) — b cot(7ibx) , 
® b (x + 1/6) = @ b (x) - — cot(nx/b). 



(C.5) 
(C.6) 



On the other hand, as follows from (JA.6J) . the function — - ^logSfe(x) satisfies exactly the same 
properties. This allows to conclude (for non-rational b) that 



@ b {x) = -- — logS 6 (z) + c , 
7T ax 



(C.7) 



where the constant c can in principle be found, but in any case our results are independent of it. 
From the property ()C.7|) . it is clear that ()C.3j) can be represented as 



sinh(7rs6) sin [f (2(3 - £)] ds 



so that the undetermined constant Cq does not enter the result. The argument of the logarithm 
coincides in fact with the factor C% i n (JDIEJ)- Therefore, the full three-point function can be 
rewritten as 



C(b,/3,(3\S,S,S') = 7T/i 7 (6 2 )6 



> 2\t2-2f> 2 



1 ^(i-^ r (2/36 - 6 2 ) r - £) ac : 



(i) 



r(i -6 2 ) 



7r sinh(7rs6) 



S 6 (2/9)-^-. (C.9) 
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Figure 8: Shift of the integration contour. 



Using the fact that 



S b (2(3) 



2tt 



0(1 



r 6 (2/5 - Q) 



r (2/36 - V) r (f - i - i) r 6 (Q - 2/3) 

and rewriting the derivative in terms of the boundary cosmological constant, one finds 

(Q-2f3)/2b 



(CIO) 



C(b,i3,{3\s,s,s') 



_d_ 



2(2(3 -Q) Tx^(b 2 )b 



u2\l,2~2b 2 



T b (2(3-Q)T- 1 (Q~2(3) 



S b ((3 + S b (P - i SJ f) Sb (P + i^) Sb (P - i^f) 



(c.11; 



Thus, the integral with respect to — fi B can be trivially evaluated and one obtains a result propor- 
tional to the boundary reflection amplitude ()2.23|) . as shown in ([2.24)1 . Note that the proportionality 
factor is very similar to the one appearing in the relation between the corresponding bulk quantities 

(HQ). 

Let us now return to the discussion of the integral in (jC.ld|) . Although, as we discussed 
above, our procedure to compute it was not rigorous, the result that we get for the two-point 
function is perfectly sensible. Therefore, one may wonder if there is an alternative way to justify 
the computation. For this we will consider C3 as an analytic function of its parameters, the 
boundary variables s, s' and momentum (3, and derive some difference equations satisfied by the 
integral. Such equations can be obtained by shifting the parameters in various ways. In principle, 
at least for generic b, these "shift equations" should fix the solution completely. Then we prove that 
the equations we derived are satisfied by the answer coming from our computation by residues. 

This strategy will turn out to be extremely useful in other cases as well: in the next part of this 
appendix, we will use it in order to perform the computation of the two-point function, starting 
from the three-point function, in the case in which fa = b or (3^ = b. In appendix ID. 11 the same 
strategy will allow us to obtain the integration rules for the three-point function that we use in the 
main text. 
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To find some shift equations for the integral (|C.ld|) . it is useful to rewrite it as 



,(i) 7r [ j± cot(irbt) — cot[irb(is + 1)] 



sinh(vr6s) C\ } = - / dt — ^ r"7i < rr- ( c - 12 ) 

A J sin[f (/5 + z^ + t)]sin[f (i-/5 + ^ + t)] 

— ICO 

Now let us consider, for instance, the shift is — > is + 2/6. It has an effect only on the arguments of 
the sine functions in the denominator. However, the change of the arguments can be compensated 
by shifting the integration variable t — > t — 1/6. One could think that this brings us back to the 
initial integral, but this is not precisely the case. The difference comes from the definition of the 
contour of integration described in section 12.21 Despite our shifts not being seen in the cotangent 
functions in the numerator, it is easy to see that they modify the integration contour in such a way 
that one pole of cot(7r6t) passes the contour from the right to the left and one pole of cot[iib(is + t)] 
from the left to the right. We show this schematically in fig. |H1 Therefore, in order to revert to the 
original integral, one has to add (subtract) the contribution of these poles. As a result, one arrives 
at the following equation satisfied by the integral (jC.12j) 



b ° a — l 



sinh(7r6s) C 



(i)' 

3 



Try 



+ — rTi rrji r 7 ^ttt (C.13) 

7r / 1 a \ ■ s+s' \ _.♦ hr / a i ; s-\-s' \ II \ / 



26 V si 4f(' 3 + i ^ £ )] sin [f(i-/ 3 +^)] Hf^-^^Hfi^^)] 

2bsi n [f (2/3-£)J V 



-cot -cot U 0+i*±M +cot f 



The shift is — >• is + 26 works in a slightly different way. It does not influence the denominator 
and changes only the argument of the second cotangent. Because of this, it is convenient to consider 
(jC12|) as a sum of two integrals and shift the integration variable, t — > t — 26, only in the second one. 
Then again the resulting integrals differ from the initial ones only by the contours of integration 
which pick up additional poles coming from the sine functions in the denominator. Considerations 
similar to the previous ones then lead to the following equation 



(e- 2ib9s - 1) sinh(vr6s) C 



"3 



[26 



(cot \irb(b-f3+i^)] -cot Inbh+i^-)] -cot \Trb(f3+i^) ] +cot Lbfb-p+i^)] ). (C. 14) 



In the same way, one can find other difference equations which involve shifts of s' and (3. It 
is easy to check that ()C.13|) . ()C.14|) . as well as these other equations are satisfied by the function 
given in (|C.8|) . Thus, the study of shift equations yields the same result that we reach by summing 
residues. This completes the analysis of the first of the three ways to derive the two-point function. 

C.2 The cases Bg S2 = B s b s and Bg S3 = B'l s 

Now we turn to the computation of the two-point function from the three-point correlator in the 
two remaining cases. The first one is obtained by taking fl 2 = b, Pi = = [3 and S2 = S3 = s, 
si = s', whereas the second case corresponds to the choice (3% = b, (3i = fii = (3 and si = S3 = s, 
S2 = s'. However, one can show that the expressions for the three-point function arising after these 
substitutions are related by a simple change of variables 



s-s' 



s'^-s' } t^[3 + t— Q-t, (C.15) 
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where we also showed the change of the integration variable that is needed. Since the final result, 
which is expected to coincide with ([Clip , is an even function of the boundary parameters, this fact 
ensures that the two possibilities we consider here will give the same two-point function. Therefore, 
we concentrate only on the second case, fis = b. 

With this choice of the parameters, the three-point function reduces to: 



C(p,p,b\s,s',s) = ir^(b 2 )b 



i2\t2-2& 2 



^ b M3) MS) M3) 
u l u 2 u 3 5 



(3) 
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,(3) 



s b (p+i*=*) s b {(3-i^y 

ioo 

S b (\+b-f3+i2*f+i)S b U+b-f3- 
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-Til 



S b (l+b+is'+t)s b (^+b+t)s b (l+b-0+i^+t)s b ^+2b-/3+i^+t) 

dt S b (f3+i^+t)s b U+b-/3+i^f+t) 



[irb(b-f)-i^+t)] sin[f (l-f3-i^ + t)] S b (l+b+t)S b (l+b+is>+t) 



(C.16d) 



(3) 

By shifting the parameters in Q in various ways, one can again find some shift equations 
satisfied by the integral (jC.16djl . Due to the symmetry of the integral under b — > 1/6, it is enough 
to find equations which use only shifts by multiples of b. Then the symmetry guarantees the 
existence of similar equations with b replaced by 1/6, and the final result must be symmetric as 
well. We find for instance the following three equations 



sin 



-2ibd, 



— sin [nb (b + is')] e 



-ib(d a +d s ,) 



sin [tt& ((3-b-i^f)] 
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~ sin[7rfe(2/3-fe)] 



sin[7r6(/3+i^)] b 



sin[7r6(2/3-6)] 
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(3) 



i[7rfe(/3+^)] S in[7rb(2/3-6)] S 6 (/3+^)S f) (/3+i^) ' 



S 6 2/3-6- 1) S 6 (is) 



(C.17c) 



The origin of the non-homogeneous right hand side in the third equation can be again traced back 
to the shift of the integration contour. Both shift operators e~^ d ^~ lhds ' and e lb ^ ds ~ 9s ^ transform the 
arguments of the sine functions in the integrand (jC.16d|) so that one has to pick the residues at 
the poles situated at t — i^f- + (3 — 2b and t = i^f- + (3 — b — \ in order to revert to the original 
contour of integration. The sum of the contributions of these four residues gives the right hand 
side of jniZc|) . 



One can show that the following expression for C 
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4iS b (2[3-b-l) S b {is) d , S b (l + b-[3 + i^)S b (l + b-[3 + i^) 



S b ((3 + i^)S b ((3 + i^) ds 



log 
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(C.18) 
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satisfies all of the above shift equations. When combined with (|C.16b|) and (jC16c|) . (|C.18|) re- 
produces precisely the expression (jC.9|) and thus leads to the same two-point function that was 
obtained in the first part of this appendix. 

We conclude that all three possibilities to get the two-point function starting from the three- 
point correlator give the same result ()2.24j) . which differs from the boundary reflection amplitude 
by a simple momentum dependent factor. This provides a non-trivial check of the cyclic symmetry 
of the representation ()2.20|) for the three-point function and gives strong evidences in favor of its 
correctness. 



C.3 One-point functions 

We can also derive the one-point function of the cosmological constant operator starting from the 
result ()2.24|) . Putting (3 = b and Si = s 2 = s, one obtains 



b 



r b (\-b) s b (is + b) s b (b) 



b v p/v J} r(l -b 2 ) sinh(7rs6) v ' 

The first integral with respect to —fi B gives the boundary one-point function of B^ s given in (|2.25|) . 
which we denote by W(s) since it corresponds to the resolvent in the matrix model formulation. 
The second integral gives the disk partition function with FZZT boundary conditions 



/ b 2 T (2- ±)T(l-b 2 ) f cosh ((b + ^ns) cosh((6-i)7rs) . 
Z(s) = - [n^(b 2 )} -Zr ^ ^ i K ) b) 1 + " » 1 . (C.20) 



b b 



As a last step one has to evaluate the derivative with respect to —[/,, keeping \i B constant. The 
result 

~fr = Vb ^ 7(fo2)] ^ r ^ T ( 1 -^ cosh ^ b ~ J) ™) = um (c - 21) 

shows that the one-point function obtained from the boundary three-point function, with the 
normalization that we have chosen, precisely coincides with the same quantity derived from the 
bulk-boundary correlator. This justifies the inclusion of the factor An into the initial expression for 
the three-point function. 



D Details on the boundary three-point function 
D.l Evaluating the integral 

In this appendix we study the integral ()5.39|) appearing in the calculation of the three-point function 
for degenerate momenta of the matter part. For convenience we rewrite it here: 
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First, following the method used in appendix we derive some shift equations satisfied by this 
integral. They can be obtained by applying various shift operators to the integrand and, if necessary, 
also by shifting the integration variable. As explained in detail in appendix one should take a 
special care of a possible change of the integration contour which can produce non-homogeneous 
contributions to the equations. In this way we find the following equations 10 



sinh - ™i) e" 
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where X(t) denotes the integrand (together with the numerical prefactor) and tji are its poles 
coming from the sine functions in the denominator 
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The corresponding residues of Z(t) were calculated in (|5.41j) . In the last equation they appear 
multiplied by the factor 



sin(Trbtji) 
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Then, by straightforward but tedious calculations, one can show that all equations are satisfied 
by the following function 
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where Aj is given in ()5.44j) or in (|5.51|) . Since the ratio of Sj, functions coincides with C^ 1 (cf. 
(|5.36c|) ). this is precisely the result ()5.46|) which we obtained in the main text by closing the contour 
of integration in (jD.ljl in the right half plane and summing the residues of the integrand only at 
the poles which are due to the sine functions. This gives convincing evidence of the correctness of 
the use of these "effective integration rules" . 

10 There are additional equations one can write besides (|D.2|) . but we think that the ones presented here are already 
non-trivial enough to illustrate the correctness of the result we are going to propose. 
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D.2 Explicit results for the three-point function 

This subsection contains explicit expressions of the three-point function for various low values 
of momenta rij. We consider all cases up to < 4 and satisfying the triangle inequality 

\ n 2 ~ n z\ <n>i < n 2 + n 3 . The results are expressed in terms of the one point functions Wi = W{ji). 
We will use the shorthand fa = /i s (rj), and we recall that the expression of the coefficients a n is 
given in ()5.32|) . 
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C(Pi,p 1 ,p 2 \ri,T 2 ,T 3 ) 
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sin 2 ^ 



SU1 ftW? - 2 cos ff ^-f- ( #| ( Wl + # 2 + 4 cos 2 £WWi ) - W^Wf 



snip 



snip 



sin 2 ^ 



166 2 „ 
— 2~A* 

7T 



i / sin — 

^ cos 2 £ (1 + 4 cos 2 |f) - 2 cos |f (1 + 4 cos |f ) W/ 3 - - T -£-W 3 W* 
^ sin b 2 



+ 



sin ff sin |f 
sm 2 p 



^ 2 2 (W3 + ^l) -2 cos 



o sin If 

2-7T 0^ 




(W 2 (W 3 2 + W^ 2 + W 3 Wi) - W£Wi 



+ l^!^ cos 4. cos 2| f( ^ i _ # ; 

+ /i 2 W? - Wl - —J^WsW^Wi - Wl) + 2 cos |f — J£Lw£W*(W 3 - W x ) 
\ snip sm ? 

' 1 t sin — sin — 
■fiv cos 2 £ (1 + 4cos 2 |f) (Wl - Wl) -g 2 - 62 #3^1(^3 - ^ 



166 : 



7T Z 



sin 2 p 



+ 



10246 4 



/i* cos 4 |f cos 2 |f(iy 3 -^ 



7T J 




(D.9) 



All of the above results satisfy the required properties of being real and cyclically symmetric in the 
momenta. In particular, cyclic symmetry is evident for n 2 , n 3 ) = (0, 0, 0) and (1, 1, 1), while we 
checked it for the other cases explicitly (notice that in order to check the symmetry one also has 
to interchange the chiralities of the operators and thus to take into account the fact that different 
renormalization factors must be used). 



D.3 d from C once again 

Our purpose here is to do a consistency check, namely to verify that the c = 1 limit of the three- 
point function produces the same two-point function which was derived in section 15.41 In fact, 
we can do better and check that one function follows from the other before taking the limit, but 
after the momenta degenerate in the matter sector were substituted and all special functions and 
integrals were eliminated. Then the agreement of the limiting expressions trivially follows. 

In other words, we are going to show that if one takes the expression (|5.52|) . puts (n\, n 2 , n 3 ) = 
(0, n, n), (Yi, t 2 ,t 3 ) = (r, r, r') and integrates with respect to — /t B (r), one obtains the result (|5.31j) . 
Of course, one can choose two other possibilities to get the two-point correlator similar to the ones 
in appendix IC. 21 However, they require more complicated calculations, and we proceed here only 
with the first choice. 

First, one can show that, after setting n\ = 0, n 2 = n 3 = n, (|5.52j) becomes 

{II (W{T 2 )-W{T 3 ~2mk)) n (W(n)-W(73-27rife)) ") 
— — — — \, (D.10) 
^12 1^23 1^12 1^13 J 

where the coefficient a n is defined in ()5.32j) . Then, taking the limit T\ = t 2 = r, t 3 = r', one obtains 
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n (W'(r)-W , (r'-27rtJfc)) d T f[ (w(t)-W(t' -2nik)) 

C(b,f3 n ,f3 n \T,T,r') = I k — h — ; } . (D.ll) 

The integral is easily evaluated producing 

„ n 

- J dft B (r)C(b,(3 n J n \T,T,r') = -^- J] (#(r) - W(r' - 2mk)) . (D.12) 

k=—n 

This result coincides with the two-point function d(f3 n \r, r') from (j5.31j) . This constitutes a check 
of the correctness of our intermediate calculations in the derivation of the three-point function. 
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